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Abstract 

Using self-avoiding walk model on three-dimensional Sierpinski fractals 
(3d SF) we have studied critical properties of self-interacting linear polymers 
in porous environment, via exact real-space renormalization group (RG) 
method. We have found that RG equations for 3d SF with base 6 = 4 are 
much more complicated than for the previously studied 6 = 2 and 6 = 3 
3d SFs. Numerical analysis of these equations shows that for all considered 
cases there are three fixed points, corresponding to the high-temperature 
extended polymer state, collapse transition, and the low-temperature state, 
which is compact or semi-compact, depending on the value of the fractal 
base 6. We discuss the reasons for such different low-temperature behavior, 
as well as the possibility of establishing the RG equations beyond 6 = 4. 

Introduction 

At low temperatures and in poor solvents linear polymers are in compact 
globule state, while at high temperatures and in good solvents they have 
extended coil configuration [1]. Self-interacting self-avoiding walk (SISAW) 
on a lattice is a good model for linear polymer in solvent. If — e < is 
the attractive energy per pair of nearest-neighbors monomers (not directly 
connected) and w =exp(e/fcsT) the corresponding Boltzman factor, then 
for iV-step SISAW and N 1 the following behavior of the mean-squared 
end-to-end distance (R%) is expected: 

(R%) ~ . (1) 

There is a critical value of interaction w = w c , such that for any w < 
w c critical exponent v = vsaw, which corresponds to extended coil phase, 
whereas for any w > w c exponent v has smaller value, equal to 1/d for 
d-dimensional homogeneous lattices. The value w = w c corresponds to the 



1 



collapse transition, with v = u c , between the low- and high-temperature v 
values: 1/d < v c < i>saw- 

On non-homogeneous lattices collapse transition also occurs, with the 
different values of the exponent v. Three-dimensional Sierpinski fractals (3d 
SFs) are nonhomogeneous lattices where exact real space renormalization 
group (RG) method can be applied for studying asymptotic properties of 
SISAW [2,3,4]. 

Renormalization group scheme 

Each member of the 3d SF family is labelled by b = 2, 3, 4, • • • and can 
be constructed recursively starting with a tetrahedron of base b, containing 
b(b + l)(b + 2)/6 unit tetrahedrons - so called generator Gi(b) (see fig.l). 




G\<4| 



Figure 1: First two stages of the construction of 3d SF with 6 = 4 

The subsequent fractal stages are constructed self-similarly, by replacing 
each unit tetrahedron of the initial generator by a new generator. To obtain 
the Zth-stage fractal lattice Gi(b), which we shall call the Zth order generator, 
this process of construction has to be repeated I — 1 times, and the complete 
fractal is obtained in the limit I — > oo. Fractal dimension df of 3d SF is 
equal to 

d f = ln[(6 + 2)(6 + l)b/6]/ln&. (2) 

In this case RG parameters are restricted generating functions, which 
are for the Gi(b) equal to 

A® (x, w) = A N , P x N w p , (x, w) = Y, B N)P x N w p , (3) 

N,P N,P 

where x is weight of any step of a SISAW (fugacity), and -Ajv,p (Bn,p) is 
the number of A-step SISAWs with P nearest-neighbor contacts, traversing 
the Gi(b) once (twice) - see fig. 2. For elementary tetrahedrons Go (b) one 
has 

A(°\x,w) = x + 2x 2 w + 2x 3 w 3 , B(°\x,w) = x 2 w 4 . (4) 
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For any b RG equations have the form 

n,m 

5 (m) = £#n,m (AV) n (BU) m , (5) 

n,m 

where A n , m (B n ,m) is the number of SISAW configurations consisting of n 
branches of "A®" — type and m branches of "!?(')" — type, traversing G^+i(b) 
once (twice). 




Figure 2: Examples of A and 5 type of SISAW on Gi(2) 

Analysis of the 6 = 4 case 

As 6 increases the number of possible SISAW configurations quickly 
grows, which is the reason why we don't quote the corresponding extremely 
cumbersome RG equations [4] for the studied 6 = 4 case here. Numerical 
analysis of these equations shows that there are three physically relevant 
fixed points. 

(1) The high-temperature fixed point (A* = 0.289912,5* = 0.0122393...) 
is reached for any w < w c = 5.028... and corresponding critical fugacity 
x = x*(w). Linearizing RG equations around this fixed point one gets only 
one relevant eigenvalue Ai = 8.6924..., and the critical exponent v is equal 
to u SAW = In 6/ In Ai = 0.641076.... 

(2) The low-temperature fixed point (A* = 0,B* = 0.356883... = 22" 1 / 3 ) 
is reached when starting with w > w c . There is one relevant eigenvalue 
Ai = 16, and the critical exponent is v = 1/2. 

(3) The tricritical fixed point (A* = 0.19292..., B* = 0.3388...) corresponds 
to the collapse transition, when w = w c . In this case there are two relevant 
eigenvalues: Ai = 15.4194... and A2 = 5.52427.... The critical exponent v is 
equal to v c = ln6/lnAi = 0.506756.... 

Discussion 

Values of v for SISAW on all so far studied 3d SF [2-4] are depicted 
on fig. 3. One can see that for any fractal u has its largest value in the 
high-temperature phase, and that it decreases as b grows, approaching the 
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corresponding Euclidean value Mjaw ~ 0.588 [5]. Respective v values are 
smaller at the collapse transition, also decreasing with 6, being close to the 
collapse Euclidean value = 0.5 [1]. Low-temperature behavior is more 
interesting. In contrast to the 6 = 2 case for which v = 1/df (meaning that 
polymer densely fills the underlying lattice, as on homogeneous lattices), 
critical exponent v for 6 = 3 and 6 = 4 is larger than 1/df, indicating that 
polymer is in semi-compact, rather than compact phase. This is due to the 
topological frustration, which is not the same for 6 = 3 and 6 = 4, being also 
the reason for strange behavior of the low-temperature RG fixed point: B* is 
equal to 22 -1 / 3 for both 6 = 2 and 6 = 4, whereas for 6 = 3 it is B* = oo (of 
course, A* = for any 6). It would be interesting to find out what happens 
beyond 6 = 4, but unfortunately the number of SISAW configurations grows 
extremely fast with 6. Estimated computer time required for enumeration 
and sorting SISAWs for 6 = 5 is not accessible, so one should search other 
methods for calculating v. 
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Figure 3: Critical exponent v for 3d SFs with 6 = 2,3, and 4 
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